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ABSTRACT 

In this work we continue a line of inquiry begun in Kanner et al. which detailed a strategy for utilizing 
telescopes with narrow fields of view, such as the Swift X-ray Telescope (XRT), to localize gravity wave (GW) 
triggers from LIGO/Virgo. If one considers the brightest galaxies that produce ~50% of the light, then the 
number of galaxies inside typical GW error boxes will be several tens. We have found that this result applies 
both in the early years of Advanced LIGO when the range is small and the error boxes large, and in the later 
years when the error boxes will be small and the range large. This strategy has the beneficial property of 
reducing the number of telescope pointings by a factor 10 to 100 compared with tiling the entire error box. 
Additional galaxy count reduction will come from a GW rapid distance estimate which will restrict the radial 
slice in search volume. Combining the bright galaxy strategy with a convolution based on anticipated GW 
localizations, we find that the searches can be restricted to about 18 ± 5 galaxies for 2015, about 23 ± 4 for 
2017, and about 11 ± 2 for 2020. This assumes a distance localization at or near the putative NS-NS merger 
range for each target year, and these totals are integrated out to the range. Integrating out to the horizon would 
roughly double the totals. For nearer localizations the totals would decrease. The galaxy strategy we present in 
this work will enable numerous sensitive optical and X-ray telescopes with small fields of view to participate 
meaningfully in searches wherein the prospects for rapidly fading afterglow place a premium on a fast response 
time. 

Subject headings: galaxies: statistics - gamma-ray burst: general - gravitational waves - X-rays: general 


1. INTRODUCTION 

The advent of Advanced LIGO (aLIGO - Aasi et al. 2015) 
and Advanced Virgo (AdV - Acernese et al. 2015) heralds the 
dawn of a new age of discovery in which gravitational wave 
(GW) detections will supplement traditional electromagnetic 
(EM) detections, such as those by Swift (Gehrels et al. 2004). 
We refer to aLIGO/AdV in combination as LVC. These GW 
observatories will begin operating soon at a fraction of design 
capability, and within a few years should be able to detect 
neutron star-neutron star (NS-NS) mergers out to ~450 Mpc 
and black hole-neutron star (BH-NS) mergers out to ^900 
Mpc (Aasi et al. 2013). These distances are referred to as the 
“horizon” for these events. 

It is important to distinguish between “horizon” and 
“range”. The “inspiral range”, the most commonly cited 
figure-of-merit as regards LVC sensitivity, is defined as the 
radius of a sphere whose volume equals the sensitivity vol¬ 
ume within which a signal-to-noise > 8 detection is achieved 
for a 1.4— 1.4M0 NS-NS merger, averaged over all sky loca¬ 
tions and binary inclinations. The antenna projection func¬ 
tions and associated averagings were discussed by Finn & 
Chernoff (1993) and Finn (1996). The inspiral range is not 
a hard upper limit, as face-on binary orbit mergers produce 
stronger signals (Dalai et al. 2006; Maggiore 2007; Schutz 
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2011; Nissanke et al. 2013). The maximum theoretical de¬ 
tection distance, the horizon, is 2.26 times the range (Finn & 
Chernoff 1993; see also Abadie et al. 2010). A simple calcu¬ 
lation shows that in a general population with binaries of ran¬ 
dom inclinations and positions over the whole sky, the frac¬ 
tion of aLIGO detections one expects to pick up from beyond 
the nominal inspiral range is about half of the total (Nissanke 
et al. 2010, 2013; Singer 2015). 

High signal-to-noise (SNR) detections will permit impor¬ 
tant physical parameters to be measured that are not easily 
accessible through traditional means, for instance the masses 
and spins of the merger components and the luminosity dis¬ 
tance Dl to the merger. A joint EM/GW detection would 
provide powerful constraints on the merger physics. Various 
groups have begun to examine prospects for finding an EM 
counterpart to a GW trigger (Nuttall & Sutton 2010; Hanna, 
Mandel, & Vousden 2014; Ghosh & Nelemans 2014; Fan, 
Messenger, & Heng 2015). Evans et al. (2015) consider 
strategies for X-ray observations using Swift/XRT. 

In Section 2 we present an overview of short GRBs, includ¬ 
ing putative GW rates for NS-NS mergers, and a discussion 
of “kilonova” emission accompanying the decay of radionu¬ 
clides. In Section 3 we describe a galaxy catalog - the Census 
of the Local Universe - and quantify its completeness out to 
distances relevant for aLIGO. In Section 4 we estimate galaxy 
totals concomitant with putative aLIGO error volumes for the 
next 5 yr, if one restricts attention to bright galaxies. In Sec¬ 
tion 5 we discuss potential EM observations in the optical. 
X-ray, and radio which might follow GW detection, and in 
Section 6 we sum up. 

2. ALIGO-ADV (LVC) OBSERVATIONS 
2.1. GW radiation 

Astrophysical sources that are powerful in EM detectors 
are always weak in GW if the underlying physical process 
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is nearly spherically symmetric. Supemovae are among the 
brightest EM transients but might produce peak GW strain 
amplitudes of \h\D < 10 cm, where h is the dimensionless 
strain and D the distance to source, and so be observable in 
aLIGO only to a few kpc (Ott et al. 2013, see their Fig. 14). 

In fact the lowest order contributions to GW radiation arise 
from a changing mass quadrupole moment. As such, binary 
mergers of compact objects give by far the strongest signals, 
yielding peak GW strain amplitudes \h\D ~1 km. Belczyn- 
ski, Kalogera, & Bulik (2002) utilize a population synthesis 
code to delineate a range of rates for different types of binary 
mergers. They calibrate NS-NS merger rates using known 
galactic binary pulsars. For aFIGO they infer 1 — 400 yr“^ 
for NS-NS, 9 - 400 yr'^ for BH-NS, and 0 - 8000 yr'^ for 
BH-BH. Abadie et al. (2010) incorporate the results of the 
Belczynski et al study and estimate an FVC NS-NS merger 
rate ^40 yr“^, with a range ~0.4—400 yr“^. A conservative 
lower limit to the NS-NS rate for a putative aFIGO inspiral 
range of 200 Mpc is ~3 yr“^ (Phinney 1991; for a more up-to- 
date study, see O’Shaughnessy & Kim 2010). More recently 
Dominik et al. (2013) calculate local (i.e., 0 ~ 0) rates of 
~10^ Gpc“^ yr“^ for NS-NS mergers and ~10 Gpc“^ yr“^ 
for BH-NS mergers. Given a putative FVC sensitivity volume 
~ 0.03 Gpc^ relevant for a 200 Mpc range, this translates in 
FVC detection rates of ~ 3 yr~^ and ^ 0.3 yr“^, respectively, 
similar to Phinney (1991). As indicated earlier, taking into ac¬ 
count the more realistic totals achieved in mapping from range 
to horizon would roughly double these estimates. In this work 
we restrict our attention to NS-NS mergers. 


2.2. GRBs and beamed EM radiation 

GRBs come in two flavors: long (> 2 s) and short (< 2 s) 
(Kouveliotou et al. 1993). Fong GRBs (IGRBs) are thought 
to arise from the explosion of a massive star - nearby IGRBs 
have associated supernovae. More interesting for aFIGO are 
short GRBs (sGRBs) which are thought to be due to NS-NS 
mergers for which we are placed along the binary rotation axis 
(Eichler et al. 1989; Narayan, Paczyhski, & Piran 1992). The 
fact that IGRBs are mostly spherically or axially symmetric 
makes them less interesting than sGRBs vis a vis GW detec¬ 
tions. 

The observed redshift range is from about 0.2 to 2 for short 
GRBs (sGRBs), with a mean of about 0.4. The XRT com¬ 
monly observes sGRBs up to 2 ; ~ 0.5 but in 10.5 yr of opera¬ 
tion Swift has not seen one with a measured redshift z < 0.1, 
i.e., the NS-NS merger aFIGO target GW horizon ~400 Mpc, 
suggesting that such events may be quite rare. 

Estimates for jet beaming are ^ 5° for IGRBs and 
0j ~ 5 — 15° for sGRBs (Burrows et al. 2006, Grupe et al. 
2006, Fong et al. 2012). Beaming angles for sGRBs are still 
highly uncertain. The beaming factors /{, = 1 — cos6*j ~ ^f/‘^ 
are roughly 1/300 for IGRBs and 1/30 for sGRBs. Based on 
the observed rate of sGRBs by Swift, Coward et al. (2012) 
estimate a FVC detection rate of ~3 — 30 yr“^ for Oj ~ 15°. 
Chen & Holz (2013) claim 3 - 7 yr"! for GRB GWh-EM 
detections. Kelley et al. (2013) estimate the rate of Swift 
or Fermi observations joint with FVC detections to be ~0.07 
yr“^. Siellez, Boer, & Gendre (2014) consider current and fu¬ 
ture high energy missions and estimate a rate of simultaneous 
GWh-EM detections of ~0.1—4 yr“^ in the FVC era. Wander- 
man & Piran (2015) estimate a co-detection rate YMC+Fermi 
of 0.1 - 1 yr-i and INC+Swift of 0.02 - 0.14 yr"!. 


2.3. Kilonova emission 

Various groups have explored the supernova-like tran¬ 
sient powered by radioactive decay of the spray of material 
~10“® — 1O“^M0 ejected from the NS (Eichler et al. 1989; 
Fi & Paczynski 1998; Metzger et al. 2010; Metzger & Berger 

2012) . The resultant “kilonova” (dimmer than a supernova 
and brighter than a nova) would produce relatively isotropic 
optical/NIR emission after a NS-NS/BH-NS merger. While 
SNIa light curves are powered primarily by decay of ®®Ni, 
the ejecta from a disrupted NS is neutron rich and yields little 
Ni. Much heavier radioactive elements form via rapid neu¬ 
tron capture (r—process) nucleosynthesis following the de¬ 
compression of the ejecta from nuclear densities. These newly 
synthesized elements undergo nuclear fission, a and /3 decays 
on much longer time-scales. The resulting energy release can 
power detectable thermal emission once the ejecta expands 
sufficiently that photons can escape. Recent general relativis¬ 
tic NS-NS merger simulations (Bauswein, Goriely, & Janka 

2013) indicate that a small fraction of the ejecta, ~1O“^M0, 
or a few percent, expands rapidly enough for most neutrons 
to escape capture. The /3-decay of these free neutrons in 
the outermost ejecta powers a precursor to the main kilonova 
emission peaking on a timescale of hours after NS-NS merger 
(Metzger et al. 2015). For D ~ 200 Mpc this emission peaks 
in the 17—band (~^365 nm) at mu ~ 22 (Metzger et al. 2015). 

Kasen et al. (2013) argue that the opacity of the expanding 
r—process material is dominated by bound-bound transitions 
from those ions with the most complex valence electron struc¬ 
ture, i.e., the lanthanides (Tanaka & Hotokezaka 2013; Gross- 
man et al. 2014; Wanajo et al. 2014). They compute atomic 
structure models for a few representative ions in order to cal¬ 
culate the radiative transition rates for tens of millions of lines, 
and find that resulting r—process opacities are orders of mag¬ 
nitude larger than that of ordinary (e.g., iron-rich) supernova 
ejecta. The resultant light curves should be longer, dimmer, 
and redder than previously thought. The spectra have broad 
absorption features and peak in the IR (~1 /rm). Kasen et al. 
(2015) combine two-dimensional hydrodynamical disk mod¬ 
els with wavelength-dependent radiative transfer calculations 
to generate model light curves and spectra. They discern two 
components to the kilonova light curve, a blue optical tran¬ 
sient (~2 d) arising from the outer lanthanide-free ejecta and 
an IR transient (^10 d) coming from the inner, lanthanide 
line-blanketed region. There had been an earlier suggestion 
of these two components in work by Barnes & Kasen (2013) 
using a less sophisticated model. 

Specific predictions for kilonova light curves are depen¬ 
dent on uncertainties such as the type of ejecta (dynamical 
and/or disk outflows), ejecta masses, and velocities. Recent 
time dependent calculations find ejecta masses in the range 
~10“^ - 1O“^M0 and velocities ~0.1-0.3c (e.g., Foucart et 
al. 2011; East & Pretorius 2012; Piran et al. 2013; Kyu- 
toku et al. 2013; Hotokezaka et al. 2013; Fernandez & 
Metzger 2013; Sekiguchi et al. 2015; Foucart et al. 2015; 
Fernandez et al. 2015). Assuming iron-rich supernova ejecta, 
Metzger & Berger (2012) predict peak optical luminosities 
~10^^-10^^ erg s“^ and concomitant values —14 to 
-17. For Meject = 1O“^M0 and Uejecta = 0.1c, Barnes 
& Kasen (2013) calculate a peak absolute magnitude in the 
near-IR (A ~ 1.7/im) of Mh = —15.5, in good agreement 
with Tanaka & Hotokezaka (2013), Tanaka et al. (2014), and 
Grossman et al. (2014). 

GRB 130603B might be the first detected kilonova (Tan- 



3 


vir et al. 2013; Berger, Fong, & Chomock 2013). It was a 
short GRB at z = 0.356 with a duration ^0.2 s in the BAT. 
If correct, it would confirm that compact-object mergers are 
the progenitors of short GRBs and also the sites of significant 
production of r—process elements. 

In addition to the optical/near-IR kilonova emission, one 
also expects a characteristic signal in the radio (Nakar & Pi- 
ran 2011; Hotokezaka & Piran 2015) as the ejecta first in¬ 
teract internally and then externally with the ISM. The latter 
interaction gives rise to a blast wave with concomitant en¬ 
hancement of magnetic fields and electron acceleration, lead¬ 
ing to synchrotron radiation and radio emission. Three tem¬ 
poral components to the radio band have been considered 
and studied, (i) early-time anisotropic emission along the rel¬ 
ativistic jet axis associated with the ultra-relativistic ejecta, 
(ii) mildly relativistic, quasi-isotropic emission accompany¬ 
ing cocoon-breakout, leading to potential radio flares for off- 
axis observers, and (iii) late-time sub-relativistic dynamical 
ejecta producing radio flares on time scales of years. The lat¬ 
ter emission should be nearly isotropic, and provide standard 
calorimetry on the global energetics of the initial explosive 
event, just as has been the case for long GRBs. 

3. CENSUS OF THE LOCAL UNIVERSE CATALOG 

Any given galaxy catalog is generally not optimal for GW 
follow-up studies. Consider two extremes: The Two Micron 
All Sky Survey (2MASS) survey (Skrutskie et al. 2006) has 
good coverage in both the northern and southern skies, but 
does not go very deep (Huchra et al. 2012). The Millenium 
Galaxy Catalog, comprising spectroscopic redshifts of galax¬ 
ies from 2dF or SDSS, is deep, but covers only a small slice 
along the celestial equator (Driver et al. 2005). An attempt 
to overcome these limitations led to the Gravitational Wave 
Galaxy Catalogue = GWGC (White, Daw, & Dhillon 2011). 
Its only limitation for our current study is that it does not 
extend beyond 100 Mpc. One of us (Kasliwal et al. 2015, 
in preparation) has amassed a catalog based on the union of 
several existing catalogs - the Census of the Local Universe 
(CLU) - which is suitable for GWh-EM follow-up studies. As 
we shall show, for bright galaxies the CLU is complete out to 
the anticipated aLIGO GW inspiral range for NS-NS mergers 
up to 2020. 

In order to show completeness, we must adopt a model for 
galaxy number density in the local universe. The Schechter 
luminosity function (Schechter 1976) provides a useful de¬ 
scription of the space density of galaxies as a function of their 
luminosity, pgi,\{x)dx = (j>*x°‘e~^dx, where x = LjL* and 
L* is a characteristic galaxy luminosity where the power-law 
form of the function truncates. It has proven to be applicable 
over up to 10 magnitudes in deep surveys (e.g.. Bonne et al. 
2015). The CLU catalog is amassed from many different sur¬ 
veys. One of the CLU data columns, htc, consists of apparent 
B—magnitudes ms for entries where they are available, and 
pseudo-ms values for sources from other bands, for instance 
2MASS. Hence for this work (j)*, L, L* —>■ <j)B*, Lb, Lb*- 
Physically, Lg represents the turn-over in the distribution be¬ 
tween a power-law for low x and an exponential for high x. 

We adopt the following values derived from B—band mea¬ 
surements of nearby field galaxies: (j)* = (1.6±0.3) X 
Mpc-3, a = -L07±0.07, L% = (1.2±0.1) x 10^°Lb,q, 
with a corresponding Mg = —19.7 ± 0.1 -I- Slogj^gh = 
—20.47 (e.g., Norberg et al. 2002; Liske et al. 2003; Gonzalez 
et al. 2006, and references therein). By comparison, for the 
Milky Way galaxy Mb = —20.42. This is based on a Milky 
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Fig. 1.— Sky maps of the CLU catalog. Shown are all 144214 galaxies 
{upperpanel), and only the 27559 galaxies for which L > L*^ {lowerpanel). 
The swath of incompleteness in both panels represents the galactic plane, 
which is excluded in many surveys. 

Way B—band luminosity 2.3 x 1 O^°Lb, 0 (Carroll & Ostlie 
1996), where Lb,q is the solar B—band luminosity, and an 
absolute solar B magnitude of 5.48 (Allen 1973). We take 
h = 0.7 based on the latest weighted overlap between the 
Planck results and the rest of astronomy (Ade et al. 2014). 

Integrating over luminosity gives integrated number density 
Po,gai = /“ PE<^t{x)dx. Although po.gai -> oo as Xi -> 0 for 
a < — 1, the integrated luminosity density diverges only for 
a < —2. One has 



(j)*L*x°‘^^ exp(—x)(ix 




( 1 ) 


where F is the incomplete gamma function. For a = — 1 
the total luminosity density is (j)*L*T{2 -f a) = (jfL* = 
1.9x 10 ®/iLb(©) Mpc“^. Dividing by a Milky Way B—band 
luminosity yields a density ~6 X 10 ^ MWE Mpc where 
MWE = Milky Way equivalent galaxy. Eor a = — 1, half 
of the luminosity density is contributed by galaxies with 
Xx /2 > 0.693. Eor the power law of interest in this study, 
a = —1.07, the cutoff lies at x^j^ > 0.626, or Mb 1/2 = 
— 19.97. This corresponds to ~0.66 of the Milky Way lu¬ 
minosity. To arrive at this Xif2 value we used the fact that 

x°'*^^e~^dx = 1.04559 for xi = 0 and a = —1.07, and 
half this value 1.04559/2 is achieved for xi = 0.626. 

Eigure 1 presents sky maps of the CLU catalog, showing 
all the galaxies, and also those for which a; > 1, where 
X = Lb/Lb*- The dark strips evident in the top panel in- 
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Fig. 2.— A comparison of the completeness measures of the CLU catalog {black) with GWGC {blue) and the 2MASS redshift survey {green) is given by 
showing frequency histogram distributions for 12 distance slices {solid) versus x = L/L*, using bins of width 0.1 dex in x. For CLU and GWGC x = Lq! Lq* , 
whereas for 2MASS x = . Shown also is the Schechter function cf)*x^e~^ AVdx {dotted), where the volume element AV is that for the given distance 

slice. For the 2MASS data, which are in the iC—band, a = —0.9 and Xi ^2 — 0.790. The vertical line segments in each panel indicate x = x-^j^- 
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dicate individual deep surveys which make up the CLU. Re¬ 
stricting the sample to intrinsically bright galaxies cleans it up 
considerably and reveals large scale structure. 

Figure 2 provides a measure of the completeness as gauged 
by the Schechter function. We compare the CLU catalog with 
the GWGC and 2MASS redshift survey out to 200 Mpc us¬ 
ing 12 radial slices. Within each slice we bin the data in 
X = LjL* and compare with the Schechter function weighted 
by the volume of the given slice. For CLU and GWGC 
X = Lb/Lb* whereas for 2MASS x = Lk/Lk*- Since 
GWGC ends at 100 Mpc there are no data in the more distant 
bins. For all three catalogs there is a progressive loss of fainter 
galaxies with distance relative to Schechter. For all radii the 
CLU follows the Schechter function for x > Xi/ 2 - In some 
of the panels in Figure 2 one sees spikes at a; ~ 10. These 
may be due to misidentified local objects spuriously placed at 
larger distances. 



Fig. 3.— Completeness relative to the Schechter function for x > a: 1/2 for 
CLU {black), GWGC {blue), and 2MASS {green) using the data presented in 
Figure 2. 

In this paper we propose to utilize the brighter galaxies in 
each bin, which tends to mitigate the effects of incomplete¬ 
ness which are brought about by a progressive loss of fainter 
galaxies with increasing distance. The CLU catalog is fairly 
complete relative to Schechter above Xi/ 2 - This is shown in 
Figure 3, where we compare with GWGC and 2MASS. In this 
work we will argue for including only the brighter galaxies in 
GWh-EM follow-up studies. As noted earlier, this implies cut¬ 
ting galaxies below Xi/ 2 - For the three fiducial distances of 
interest in the next section, 60 Mpc, 120 Mpc, and 180 Mpc, 
the CLU catalog is complete above this cut line at a level 
^100%, ~80%, and ~40%, respectively. Since these three 
distances are based on the NS-NS inspiral ranges for 2015, 
2017, and 2020, an additional factor which must be included 
for the third year stems from the fact that the CLU catalog 
ends at 200 Mpc, namely the fact that about half the detec¬ 
tions would be expected to arise from between the range and 
the horizon. Therefore for the third year the completeness is 
^20%. Both Figures 2 and 3 show some regions of “over¬ 


completeness” below 100 Mpc. This is probably due to local 
over-densities such as the Virgo Supercluster. 

Our choice of Xcut = a:i /2 is motivated by a trade-off be¬ 
tween including enough galaxies so that we encompass a rea¬ 
sonable fraction of putative sGRB host galaxies (e.g., Berger 

2014, see his Figure 8) on the one hand, and not having in¬ 
completeness become too great an issue on the other hand. 
Taking Xcut = X\i 2 picks out about half of the sGRB host 
galaxies shown in Berger (2014), therefore the true efficiency 
factors for the three target years would be ^50%, ~40%, and 
^10%, respectively. 

Another caveat arises from the fact that by utilizing 
B—band data we are implicitly taking B—band luminosity 
as a proxy for the compact object merger rate. Although 
this assumption is generally supported by sGRB observations 
(Berger 2014), it may in fact not be a universally good proxy 
(e.g., de Freitas Pacheco et al. 2006; O’Shaughnessy et al. 
2010; Hanna et al. 2014). 

4. GALAXY STRATEGY 

As regards identifying a candidate host galaxy for a GW 
event, a galaxy catalog is only half of the picture. Of course 
the starting point is the GW localization. Initial localizations 
in the ~2015 time frame are expected to be ~500 deg^, de¬ 
creasing to ~20 deg^ by ^2020 (Aasi et al. 2013; Singer et 
al. 2014; Berry et al. 2015). In terms of the depth, the NS- 
NS inspiral range should increase from ~60 Mpc to ~180 
Mpc over the same time frame. A projected timeline of NS- 
NS inspiral range versus date is given by Aasi et al. (2013). 
One may obtain a good estimate using simple considerations. 
In the CLU catalog there are N* = 27559 galaxies brighter 
than Mb*, and V 1/2 = 47438 brighter than Mb 1 / 2 - If one 
restricts A^i /2 based on the limited sky areas and volumes rel¬ 
evant for the three target years, i.e., 500 deg^ x 60 Mpc for 

2015, 100deg2 x 120 Mpc for 2017, and20deg2 x 180 Mpc 
for 2020, one obtains A"gai = 26.0, 36.4, and 18.8, respec¬ 
tively. As noted previously, since these distances represent 
putative NS-NS inspiral ranges, for a localization near fi we 
must roughly double the totals to take into account going from 
range to horizon if the distance error is large. 

Realistic idealizations for aLIGO localizations have been 
undertaken by several recent groups which attempt to quan¬ 
tify the “volume reductions” which might realized. Nissanke 
et al. (2010) use a Markov Chain Monte Carlo (MCMC) 
technique to calculate distance measurement errors associated 
with aLIGO localizations of astrophysical populations of NS- 
NS and BH-NS binaries, considering both isotropically ori¬ 
ented as well as beamed events. They take as their starting 
point a precise sky localization based on a coincident EM de¬ 
tection of the same GW event. They present Eisher-matrix- 
derived linear scalings for [AB^/True Db] for the two pop¬ 
ulations, assuming four GW detectors. If the EM emission 
from the NS-NS merger providing the coincident EM signal 
is isotropic, they find that, in combination with the precise 
sky position, the distance to NS-NS binaries can be measured 
with a fractional error of ^20 — 60%, with most events clus¬ 
tered near ^20 — 30%. If the EM emission from the NS- 
NS merger is beamed, with a ~25° opening angle, then the 
error on the distance is reduced by a factor ^2 and much 
of the high error tail is eliminated. BH-NS events are mea¬ 
sured more accurately: the distribution of fractional distance 
errors lies in the range ~15 — 50%, with most events clus¬ 
tered near ^15 — 25%. Nissanke, Kasliwal, & Georgieva 
(2013) carry out extensive end-to-end simulations, looking at 
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GW sky localization, distance errors and volume errors using 
NS-NS and BH-NS mergers. They compare MCMC-derived 
distance measures with three dimensional (3D) volume mea¬ 
sures. They outline optimal strategies to prepare for identify¬ 
ing EM counterparts of a GW merger. 

Singer et al. (2015) have created full 3D position recon¬ 
structions for a large population of simulated early aLIGO 
NS-NS events. They provide a simple approximation for the 
3D distance distribution and qualitatively describe the shapes 
that emerge. In the present work, as we are more interested 
in the impact of the galaxy catalog we use an even simpler 
description; we assume that the localization subtends a given 
solid angle, and is a shell between two constant radii. 



Fig. 4.— The number of CLU galaxies in a given size error box versus 
range {black). In this experiment we implicitly assume a localization near 
the inspiral range for each of the three target years. Diamonds indicate the 
adopted NS-NS inspiral ranges for the three target years of this study, /i = 
60, 120, and 180 Mpc, respectively. The three curves are representative of the 
increasing localization capability of aLIGO+VIRGO with time, AQ ~ 500 
deg2 (for 2015), 100 deg^ (for 2017), and 20 deg^ (for 2020). The CLU 
incompleteness weighting factors at r = fi are 1.0, 1/0.8 = 1.25, and 
1/0.4, respectively. The galaxy count totals indicated by the diamonds are 
^gal = 17.7± 5, 22.5 ±4, and 10.5 ± 2. We also present the corresponding 
Afgai values determined directly from the Schechter function (red) weighted 
by the relevant volume AU = (4/3)71/i^(Af2/47r) (see the end of Section 

4). 

Singer et al. (2015) calculate full 3D aLIGO reconstruc¬ 
tions based on BAYESTAR, a rapid reconstruction method 
for BNS mergers, and LALInference, a full Bayesian param¬ 
eter estimation code (Singer et al. 2014; Singer 2015; Veitch 
et al. 2015), for the near-future of aLIGO - 2015 and 2016. 
In this work we consider a longer time frame, i.e., up to the 
full achievement of aLIGO design sensitivity. In order to cal¬ 
culate galaxy sky counts in error boxes of given areas on the 
sky, we take a simplified approach compared to Nissanke et al. 
(2009), Nissanke et al. (2013), and Singer (2015). Namely, 
rather than using a realistic 3D reconstruction, for a given 


line-of-sight (LOS) we take a simple top-hat windowing func¬ 
tion. Lor our three putative aLIGO target years - 2015, 2017, 
and 2020 - we adopt p, = 60, 120, and 180 Mpc, respectively, 
as fiducial LVC NS-NS inspiral ranges. 

We consider 1000 randomly selected LOSs over the sky, 
and then search the CLU catalog to find galaxies within an 
angular separation that would place them inside an error box 
on the sky of (i) 500 deg^, (ii) 100 deg^, or (iii) 20 deg^ for 
the three cases. Although the actual sky-projected aLIGO 
localizations will be complicated, the more important fac¬ 
tor is simply the total sky area involved. Only galaxies are 
considered for which their luminosity places them above the 
50th percentile mark Xi/ 2 - The CLU galaxy count in a ra¬ 
dial bin is incremented only if |(r — p)/p| < 0.3 (Hanna et 
al. 2014), i.e., if the galaxy lies within a thick spherical shell 
with Ap/p = 0.6, where the NS-NS inspiral range p = 60, 
120, and 180 Mpc, respectively. This has the effect of es¬ 
sentially doubling the galaxy counts derived by truncating the 
integration at the range p, thereby mimicking the effect of the 
range-to-horizon mapping (Hanna et al. 2014). Thus in this 
experiment we are implicitly assuming a localization near p 
for each target year, and thus in some sense our galaxy count 
totals represent upper limits. 

Ligure 4 shows the results of this experiment. Each point 
is the mean of the 1000 individual LOS values, and the er¬ 
ror bar is the standard deviation. Two main points are worth 
noting. The most obvious is simply that larger error boxes 
yield more galaxies. The second is that counts are added only 
out to roughly the range for a given year. The CLU totals 
out to r = p (for a localization in each year near p) are 
A"gai = 17.7 ± 5, 22.5 ± 4, and 10.5 ± 2, respectively. As 
noted previously, the CLU is ^100% complete above 
60 Mpc, ^80% complete at 120 Mpc, and ~40% complete at 
180 Mpc. In addition, the NS-NS range-to-horizon mapping 
roughly doubles the totals for localizations near r = p with a 
large radial uncertainty. 

Lor consistency we may consider galaxy counts derived di¬ 
rectly from the Schechter function. The galaxy density above 
4^1/2 is (j>*x°‘e~^dx = 2.35 x 10“^ Mpc“^. If we then 

multiply this by a volume AU = (4/3)7rp^(AU/47r), where, 
for the three target years (p, AO) = (60 Mpc, 500 deg^), (120 
Mpc, 100 deg^), and (180 Mpc, 20 deg^), respectively, we 
obtain A"gai = 25.8, 41.3, and 27.8. 

5. TILING OBSERVATIONS 

The results given in Ligure 4 reveal the dramatic reduction 
in tiling requirements brought about by restricting a search 
methodology to the brighter galaxies. Were one simply to tile 
the entire 2015, 2017, and 2020 putative error boxes A = 500, 
100, and 20 deg^ using an EM detector with a field-of-view 
(LOV) of 5A = 0.1 deg^, i.e., ~19 x 19 arcmin, the number of 
tilings required would be AjbA = 5 x 10^, 10^, and 2 x 10^, 
respectively. 

Concentrating on selected bright galaxies reduces the tiling 
effort considerably. The optimal size of a “tile” is dictated 
by observations of short GRBs - in particular their locations 
within their host galaxies. Long, Berger & Lox (2010) use pre¬ 
cise HST localizations to study the cumulative distribution of 
projected physical offsets for short GRBs with sub-arcsecond 
positions and find them to lie within ~100 kpc of their host 
galaxy centers. Lor a typical distance of interest in this study, 
100 Mpc, this corresponds to ~10“^ radian or ~0.057 deg - 
a projected area on the sky of ^0.010 deg^ or ~37 arcmin^. 
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This can be covered by an EM detector with a rectangular 
field-of-view (FOV) of at least 7x7 arcmin or 49 arcmin^. For 
such detectors the number of tilings will reduce to the number 
of galaxies as shown in Figure 4, i.e., ~20, ~20, and ~10, re¬ 
spectively. This represents a reduction in requisite tilings by 
~1 — 2 orders of magnitude over a brute force methodology, 
i.e., simply covering the entire GW error box. 

Furthermore, the underlying strategy of focusing attention 
only on bright galaxies is strengthened by the observation that 
short GRBs tend to lie in the larger, brighter galaxies (in con¬ 
trast to long GRBs which lie preferentially in dwarf irregu¬ 
lars). Fong et al. (2010) compare the projected physical off¬ 
set distribution between short and long GRBs and find that, 
when normalized to the sizes of their host galaxies, the dis¬ 
tributions are indistinguishable. However, the absolute length 
scales differ by a factor of ~5 (see also Fong & Berger 2013). 
Berger (2014, see his Figure 8 ) plots all known sGRB host 
galaxy x = Lb/Lb* values; they span a range 0.1 < x < 2. 
Our adopted cut value in this study Xi /2 — 0.6 is roughly at 
the median of the observed distribution, therefore our bright 
galaxy strategy would pick up about half the sGRBs shown in 
Berger’s sample. Therefore the total effective completeness 
would be reduced by another factor ~2 below that given in 
the previous section, which only considered completeness for 
X > XYji- The trade-off against lowering our cut value so as 
to include a greater fraction of putative sGRB host galaxies, 
i.e., using for example Xi /3 or a;i/ 4 , is two-fold: (i) the num¬ 
ber of galaxies would increase rapidly and (ii) incompleteness 
would become a more severe issue. 

5.1. Optical 

Within an area of several square degrees on the sky there 
are many more optical transients than in other wavelengths, 
and therefore more opportunities for false positives (Kulka- 
rni & Kasliwal 2009; Drout et al. 2014; Tanaka et al. 2014; 
Cowperthwaite & Berger 2015; Singer et al. 2015). Most 
of this activity stems from variable sources within our own 
galaxy. If we restrict our attention to bright galaxies, false 
positives could still arise from powerful sources within the 
target galaxy such as supernovae (SNe). However, SNe are 
rare and also have longer timescales than EM counterparts to 
NS binary mergers. More generally, Nissanke et al. (2013) 
and Kasliwal & Nissanke (2014) show that false positives 
can be due to both foreground stars (e.g., flare stars and 
cataclysmic variables) and background galaxies (e.g., super¬ 
novae and AGN). Chance associations of such transients with 
a host galaxy location would be problematic to exclude based 
solely on photometry, and would probably necessitate multi¬ 
wavelength observations. An advantage with the optical (and 
near-IR) is that the emission of interest is due to the kilo- 
nova, which is quasi-isotropic, whereas X-ray emission, for 
instance, would be beamed. 

A variety of large optical telescopes that have been ac¬ 
tive in GRB follow-up have large FOVs which would make 
them amenable to GWh-EM tiling observations: GTC/OSIRIS 
(10.4m aperture - 7.8 x 8.5 arcmin FOV), Keck/FRIS (10m 
- 6 x 8 arcmin), FBT/FBC (8.4m - 23 x 23 arcmin), and 
VFTWIMOS (8.2m - 14 x 14 arcmin). These instruments 
could all cover the Fong et al. (2010) short GRB projected 
area in one tile. 

5.2. X-ray 

The transient sky is less chaotic in X-rays than in the opti¬ 
cal. Therefore individual transients stand out more. However, 


for short GRBs only the beamed events will be detected in X- 
rays, which reduces the detection chances by the reciprocal of 
the beaming factor /& of a short GRB, which is highly uncer¬ 
tain, ~ 10 — 100. The primary instrument of relevance 
for X-ray follow-up is the SwiftfXRT (Burrows et al. 2005), 
a focusing X-ray telescope with a 110 cm^ effective area, an 
18 arcsec resolution (one-half power diameter) in the 0.2 — 10 
keV band, and a field-of-view (FOV) of 23.6 x 23.6 arcmin, or 
~0.15 deg^. Thus one XRT tile would be 5A ~0.15 deg^, or 
about ten times the minimum required to tile a putative short 
GRB host galaxy. Given an XRT exposure time of ~0.1 ks 
(Kanner et al. 2012) and a comparable time to slew between 
tiles, a complete search of the three error boxes depicted in 
Figure 4 would place modest demands on Swift. The bright 
galaxy strategy for GW-EM follow-up described in this work 
results in a far less strenuous use of XRT resources than dis¬ 
cussed for example in Evans et al. (2015). 

5.3. Radio 

Explosive transients which eject ionized matter into a sur¬ 
rounding medium eventually produce radio waves via syn¬ 
chrotron radiation as the electrons interact with tangled mag¬ 
netic fields in shocked regions. Radio emission has been seen 
in GRBs, supernova remnants, colliding winds in massive bi¬ 
naries, symbiotic stars, and cataclysmic variables. Over long 
timescales the radio provides good calorimetry on the total 
energetics of the explosion due to its quasi-isotropic emission 
(Nakar & Piran 2011; Hotokezaka & Piran 2015). Thus, as 
with the optical observations of kilonovae, the radio band also 
has an advantage over X-rays, which are beamed. The radio 
sky is also relatively quiet as regards fast transients, although 
this may be at least partly due to the lower sensitivity of tran¬ 
sient radio surveys compared to optical or X-ray ones. (Frail 
et al. 2012; Mooley et al. 2013, 2015). 

Radio facilities have done extensive follow-up work on 
GRBs and would be relevant in radio tilings of GW error 
boxes insofar as having beam FWHMs that would cover the 
expected region in an F ~ L* galaxy: e.g., AMI-FA: 6 ar¬ 
cmin beam at 16 GHz and JVFA: 9 arcmin beam in C band 
(4.5—5 GHz). At least 4 short GRBs have had radio af¬ 
terglows detected, 051221a, 050724, 130603b, and 140903 
(Chandra & Frail 2012; Berger 2014). They were all quite 
faint, ^100 /iJy, but they were also considerably beyond the 
range of current interest d < 200 Mpc 01 z < 0.047. 

6. DISCUSSION AND CONCLUSION 

In this work we have argued for a strategy in which only 
the brighter galaxies are considered in GW-EM follow-ups. 
We show that the CFU catalog is fairly complete for x > 
Xi /2 out to 200 Mpc. By weighting the galaxy counts within 
projected sky areas with a simple model for the aFIGO radial 
localization, we find that only about 18 ± 5 (for 2015), 23 ± 4 
(for 2017), or 11 ± 2 (for 2020) galaxies need be considered, 
assuming error boxes of 500, 100 and 20 deg^, respectively. 
This results if one restricts attention to galaxies in the upper 
50th percentile in integrated luminosity density. Furthermore, 
there are numerous EM detectors with the property that one 
tile (i.e., one FOV) would encompass the region of interest. 
These facilities - optical. X-ray, and radio detectors - have 
carried out GRB follow-up previously and could participate in 
future tiling observations. Having one tile per galaxy reduces 
the number of requisite tiles down to the number of galaxies, 
a reduction in tiling effort by ~1 — 2 orders of magnitude. 
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The CLU efficiencies for x > Xi /2 vary from ^100% at 
NS-NS inspiral range /r = 60 Mpc to ~40% at p. = 180 
Mpc. For our third target year, 2020, the efficiency is cut an 
additional factor of two due to the fact that ~half of the detec¬ 
tions at any given time come from beyond the range, and the 
CLU is truncated beyond 200 Mpc. In addition to the CLU 
incompleteness, our adopted cut Xcut = a;i /2 lies at roughly 
the median in the observed sGRB host galaxy Lb/Lb* dis¬ 
tribution (Berger 2014), therefore we lose another factor of 
two in potential sGRB hosts. Lowering our Xcut value fur¬ 
ther would allow us to cover more of the expected sGRB host 
galaxy Lb/Lb* distribution but would exacerbate the CLU 


incompleteness issue. An additional caveat arises from the 
hostless sGRBS with optical afterglows, and the (majority) of 
sGRBs for which no afterglows have been observed. If these 
lie further from the host galaxy centers than the ~100 kpc 
maximal offset indicated by Fong & Berger (2013), then our 
strategy could miss the kilonova emission. 
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